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alloys. The method^f BozzoI^^mnte^aml^mhh^Pl^^'Rev B^s’ 4V3 uWlT ° f Cu ‘ Nl 

energetics of these systems as a function of ,„ m . Y ° 45 ’ 493 (1992) J ls used to compute the 

dictions are compared with other theoretical andTxpeTme^lults. ^ ^ C ° nCen,ratl0n ' The P"* 


Early experimental evidence obtained from work- 
function measurements showed that the surface composi- 
.on of evaporated Ni-Cu alloy films was constant, 3 
that the surface was enriched in Cu over a wide range of 

S sT P ° Sltl ? n ■ ThCSe FeSUltS were confined 

^subsequent expenmental evidence obtained in a 
broader range of temperatures, 4 and supported by hydro- 

ments a on 0 C T ZZ J b 5 e 6 nzene hydrogenation measure- 
ments on Cu-Ni films. 3 - 5 - 6 The behavior of alloy powders 

S “ K in<iica,ed made for 

lms were valid for powders as well . 7 - 9 Soon thereafter 

(AETdam^ betWee " Au 8 er " elec tron spectroscopy 

troths m C regl ° n °- 8 ~ 10 keV ’ and volumetric ad- 
sorption measurements with CO at room temperature 10 

Tr k r!' ™ by "" AES ' ,1I,erim «" S ° f Hell and c'o. 
workers. These experiments, while contradicting the 

nalysts of Takasu and Shimizu 12 that established surface 

inthe bulk d 1 e H ment Whh the min ° rity com P^ition 
Z Z f ’ , y determined the Cu enrichment of the 

it C or a "’d 6 ran 8 e of bulk compositions. These re- 
sults were later confirmed by the AES work of Kuijers 
and Ponec, and the analysis of surface segregation and 

bastofPal 1 h CC Pr °!! ieS d ° ne fr ° m 3 calcula, i° n on the 
basis of Palmberg s physical mechanism. 14 

tain 1 " 6 ™. mp ° sdion depth Profit was at this point uncer- 
tain. although most results favored Cu enrichment of 
plants below , he surface. » '* o.hers indie, ed a mZL 
decrease in Cu composition 17 or oscillations. 18 More 
precise results based on time-of-flight atoms-probe (AP- 
FIM) experiments by Ng and co-workers 19 indicated sur- 

amSrf" ° n dependence on the ^ystal face, as well 
s evidence of a monotonic decrease in Cu concentration 

In tHe SUrfaCe ’ returnin g ^ hulk composi- 

;° c i a C ^ yerS '. ° ther exper ’ me ntal techniques used 
include low-energy ,on scattering (LEIS), 20 particularly 
sensmve to the outermost atomic layer of the surface 
and x-ra> photoelectron spectroscopy (XPS) 21 The XPS 
ana^sis of Webber « al. 16 provided important evidence 
of Cu segregation in Ni-rich alloys. The XPS experi- 


ments of Wandelt and Brundle 21 also provided evidence 
for crystal face specificity. Ce 

During this same period, there was extensive theoreti- 
cal work devoted to the problem of surface segregation in 
Jj y 22 ’ lncludln 8 earJ y calculation models for binary al- 
'° yS k , the regular solut ion model, 23 calculations based on 
bulk strain energy considerations, 22 - 24 and the role of 
thermodynamte parameters of alloys. 25 - 26 Microscomc 
dectrontc theories were also applied to explain Cu enrich- 
ment, as were semiempirical methods [the embedded- 

niques^ 30< wh'^h^^ ^ 28>] a " d M ° nte Carl ° tech ‘ 
ntques, which also provided information on cluster- 

size distribution at the surface. 

More recently, Sakurai et al , 31 reported AP-FIM re- 
sults on the surface composition of Cu-Ni alloys which 
indicated Cu depletion on the surface for alloys contain- 
mg more than 84 at. % for Cu-rich alloys, 32 showing that 
o such reversal of the segregation species occurs, in ac- 
ordance with the body of previously available experi- 
mental and theoretical evidence, 1 - 30 together with ? t wo 
recent theoretical studies confirming those results 33 - 34 

fnr ,n rePOrt the , applicati °n of a semiempirical method 

Sm i ril °(^FT) C35lt 'h^t! e I? P f d by BOZZOl °’ Ferrante ’ and 
Smith (BFS), which builds on ideas underlying the 

aonliedT CrySta ‘ the0r> ' (ECT) ' 36 This method has ^been 
struefnr 2 T Cty ° f fundamen 'al Problems in alloy 

experimem r ° V H ng ’ ‘ n aH r CaSeS ’ eXCellent a greement with 
experiment. Here we focus on the study of surface 

composition and in-depth profiles of Cu-Ni alloys for 

different crystal faces, by means of a Monte Cark/calcu- 

Th" Rpc hlCh Z all ° y ener g etics are described by BFS 
The BFS method ts based on the idea that the energy of 
formation of an arbitrary alloy structure is the superposi- 

mThfano,''" “ al COn " ,bUnons ofnonequivalen, aioms 
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c, has two components: a strain energy e s that accounts 
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for the actual geometrical distribution of the atoms sur- 
rounding atom®, computed as if all its neighbors were of 
the same atomic species, and a chemical energy e .which 
takes into account the fact that some of the neighbors of 
atom i may be of a different chemical species. The cou- 
pling function g„ ensures the correct asymptotic behavior 
of the chemical energy contribution. The strain e " erg J e 
a pure defect crystal is calculated assuming y 

neighbor of atom i is of the same species X e, is then 
computed with any technique (first-pnnciples metho >d , 
semiempirical techniques, etc.). For ef we interpret the 
chemical composition as a defect of an otherwise pure 
crystal. We represent this defect by perturbing 
electronic density of the overlap region between dissimi- 
lar atoms, and locating them at equilibrium lattice ^te 
atom i. The ideas of equivalent crystal theory (ECT) are 
used to develop a procedure for the evaluation of the en- 
ergy associated with this “defect.” To free the chemica 
energy of the structural defect energy, which should be 
included only in the strain energy, we reference e, to a 
similar contribution where no such perturbation is in- 
eluded (e t C °). Finally, the coupling function g, is define 


where of 


is 


a solution of 


as gi—e » 7 . , 

e s =E‘ r [ \-(\+a?)exp(~a?)}, and where E c 1S e 

cohesive energy for atom /. We direct the reader to Ref_ 

35 where a detailed description of the calculation o 
strain and chemical energy contributions . l. 

Except for two parameters determined y ® .. 

[mental or theoVe.icai alloy properttes. the method rehes 
on pure element properties. Moreover, the pairwise cha 
acter of the interaction between dissimilar species aci i- 
tates application to multicomponent systems with no fur- 
ther experimental or theoretical input. Within the frame- 
work of BFS, the calculation of defect energies requires 
information on the atomic positions only. The input pa- 
rameters are readily available for a 

allovs The experimental input used for the Cu-Ni alloy 
t the heat of solution in the dilute hmit 
£“'[Cu(Ni)]=0.03 and £*>'[Ni(Cu)]=0.09 eX The BFS 
parameters thus obtained, together with other binary sys- 

Th; ar ;r“ in ui e d to determine the segregation 

profiles is a variant of the Metropolis Monte Carlo 
method Initially a computational cell consisting of 15 
layers of 98 atoms [(100) and (111) surfaces], or 30 layers 
of 49 atoms [(110) surface] is set up. A starting composi- 
tion is chosen, and each atom within the cell is assigned a 
speciesprobabilistically , so that the initial composition is 
uniform^ both throughout the cell and within each layer 
individually. All but the last four layers are active in 
that their compositions are allowed to vary 
computation. For reasons having to do with the deta 
of the energy calculations, the last four layers are static, 
ch^ neither composition no, the derailed dtsmbu- 
tion of the two chemical species within each layer. 

As he simulation proceeds, a pair of atoms of opposite 
siectes Ts chosen randomly from within the active region 

of the computational cell, and the total energy of the cell 

is computed. The chemicai species of the I two atom are 
reversed, and the total energy is recomputed. The 


sal is accepted according to the Metropolis criterion (it is 
accepted if it lowers the energy, or ^th probability 
e — t^E/kT if j t raises the energy by an amount ). me 
simulation is continued until the segregation pro e at 

tains a steady state. A 

Note that the procedure as described results in 
pletion of the segregating species in all but the top few 
layers. A variety of replenishment algorithms have been 
used elsewhere; in this work, however we use a computa- 
tional cell containing enough layers that the asymp 
concentration far from the surface is evident, and take 
the bulk concentration to be that value. Atomic relaxa- 
tion, which we have ignored, can afreet the result 
some cases. Previous work by the authors on ® ,ngl * 
impurity-atom segregation has shown that when there is 
a large lattice constant mismatch between the specie 
(Cu-Ni has a small lattice mismatch), the segregation en 
ergy is substantially different when computed with and 
without relaxation; the strain component of the energy 
dominant. Still, it is expected that near , the ! h ' 

where the concentration differs appreciably from J 
bulk relaxation may reduce the strain component of the 
energy, and result in a modification of the segrega ion 

profile. This issue is being investigated further. 

P Our goal is to present a set of results addressing the 
different issues covered in previous studies, concerning 
lire dependence of .he profiles op bulk 
temperature, and crystal face, as well as specific featur 
of the profiles. In agreement with experiment we 
found that there is no segregation reversal for 
concentrations as found by Sakurai et ^ or ^P^ n 
in the second layer as obtained with EAM. Moreover, 
the surface concentrations are in excellent agreement 



FIG 1 Cu(lll) surface composition vs Cu bulk composition: 
comparison between the experimental and theoretical results for 
the first and second planes. 
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with experimental values, in particular for the (111) face, 
the subject of most experimental studies. Figure 1 
displays our predictions for the (111) surface concentra- 
tion, as well as several experimental measurements and 
other theoretical approaches. Our results for the compo- 
sition of the surface plane and the first plane below the 
surface were obtained at the fixed temperature T = 800 
K. BFS predicts the formation of a solid solution for the 
bulk alloy for temperatures between 300 and 400 K. Fol- 
lowing Foiles, 28 we also computed the Warren-Cowley 
short-range order parameters 37 as a function of concen- 
tration for a wide range of temperatures. These results 
are comparable to experimental values, and clearly indi- 
cate the formation of the solid solution for temperatures 
above 400 K. Therefore, the temperature chosen for this 
calculation is well within the continuous solid solution re- 
gion of the phase diagram predicted by our model (simi- 
lar to that used in other theoretical calculations, and 
within the range of experimental measurements). Our re- 
sults seem to be particularly good for low bulk concentra- 
tions, where there is a rapid increase in the surface con- 
centration. 

Figure 2 shows the segregation profiles as a function of 
the bulk concentration and crystal face. We find a direct 
relationship between the surface concentration and the 
roughness of the surface (defined as the reciprocal of the 
surface packing fraction for hard spheres). This propor- 
tionality is also seen in the “thickness” of the surface re- 
gion, i.e., the planes near the surface where the concen- 
tration differs considerably from the bulk value. Our re- 
sults are in quantitative agreement with experiment 21 and 
other theoretical calculations. 28-30 Our results (Fig. 3) 
for the temperature dependence of the profiles agree qual- 
itatively with a previous Monte Carlo study by Donnelly 
and King. 29 The BFS values (Fig. 2) disagree quantita- 
tively with the results of Ref. 29, in part due to the fact 
that their surface layer is defined as consisting of atoms 
which are not fully coordinated. Thus, using the same 
definition for the first atomic layer in both cases, the ac- 
tual Ni composition on the first atomic layer is in good 
agreement with our predictions. Finally, Fig. 4 shows the 



FIG. 2. Cu composition depth profiles for the (111), (110), 
and (100) surface of Cu-Ni alloys. 



FIG. 3. Cu composition depth profiles as a function of tem- 
perature. 




FIG. 4. Cu composition depth profiles for different bulk con- 
centrations for the (a) (100) and (b) (110) Cu-Ni surfaces. 
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segregation profiles for two different orientations and 
several bulk concentration values. Consistent with previ- 
ous theoretical calculations, the composition approaches 
the bulk alloy within a few atomic layers, regardless of 
the bulk concentration value. On the other hand, the 
surface concentration values are consistent with those 
found experimentally for the (111) face (see Fig. 1). 

In this paper we examined several points which we be- 
lieve are of interest to the community at large: we 
demonstrated the ability of our semiempirical method to 
describe this type of phenomenon. We also provided a 
comprehensive study of the segregation profiles in terms 


of temperature, crystal face, and bulk concentration. The 
temperatures and bulk concentration values used corre- 
spond to the solid solution region of the Cu-Ni phase dia- 
gram. Our results are in agreement with several previous 
studies that indicate a monotonic profile for all tempera- 
tures and concentrations, with no evidence of segregation 
reversal. 
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